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ON SOME GENERAL SOLUTIONS OF THE SIMPLE PELL 

EQUATION 


VLADIMIR PLETSER 


Abstract. Two theorems are demonstrated giving analytical expressions of 
the fundamental solutions of the Pell equation X 2 — DY 2 = 1 found by the 
method of continued fractions for two squarefree polynomial expressions of 
radicands of Richaud-Degert type D of the form D = (/ (u)) 2 d= 2 Q: n, where 
D, n > 0, a > 0, G Z, and / ( u ) > 0,6 Z, any polynomial function of it 6 Z 
such that f (u) = 0 (mod (2“ —1 n)) or f (ti) — (2 a ~^Tlj ( mod (2 a ~ * 1 ^))- 
Keywords: Continued fraction development (11J70; 11Y65); Solutions of 
Pell equation (11D09) 


1. Introduction 


Pell equations of the general form 

(1.1) X 2 -DY 2 = N 


with X , Y. N £ Z and squarefree D > 0, £ Z, have been investigated in various forms 
since long, already by Indian and Greek mathematicians and further in the 17th 
and 18th centuries by Pell, Brouckner, Lagrange, Euler and others (see historical 
accounts in [1], [6] 1 , [7], [14]) and are treated in several classical text books (see 
e.g. [3], [9], [10], [15] and references therein). 

For N = 1, the simple Pell equation reads classically 

(1.2) X 2 - DY 2 = 1 

Beside the obvious trivial solution (X t , Y t ) = (1,0), a whole infinite branch of 
solutions exists for n > 0, £ Z given by 


(1.3) 


(A, + /DYi)” + (A, - /DFi) 


Y n = 


2 

(Ai + VDYtY - (A, - y/DYij 
2 x/D 


where (A-,, Y,) is the fundamental solution to (1.2), i.e. the smallest integer solution 
(Ai > l,Yi > 0, £ Z+) different from the trivial solution. Robertson [12] lists five 
methods to find the fundamental solution (Ai,Yi). Among these, the classical 
method introduced by Lagrange [5], based on the continued fraction expansion of 
the quadratic irrational \fD , is at the heart of several other methods (see also [8]) . 
In this paper, we demonstrate two theorems yielding analytical expressions of the 
fundamental solutions (Ai ,Y{) of the simple Pell equation for squarefree polynomial 
expressions of radicands of Richaud-Degert type D of the form D = (/ ( u )) ± 2 a n 


1 [6] lists 138 references of articles on Pell equations from 1658 to 1943. 

1 
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where D,n £ Z + , a > 0 £ Z* and f (u) any polynomial function of u £ Z with 
f (u) > 0, £ Z+ such that first, f (u) = 0 (mod (2 a ~ 1 n)) and second, f (u) = 
(2 Q_2 n) ( mod (2“ _1 n)). 


2. Solutions of the Pell equation by the continued fraction method 


Quadratic irrationals x of the form 


( 2 . 1 ) 


VD + b 

x =- 

c 


where 6, c ^ 0 and D > 0, £ Z, x £ R, and D is not a perfect square (y/D ^ Z), 
have been studied extensively (see e.g. [2], [4], [11], [13]). It is known also that a 
real number x has an eventually periodic regular continued fraction expansion if 
and only if a: is a quadratic irrational. It yields that for every squarefree D 1 there 
exists r > 0, £ Z + such that the regular continued fraction expansion of \[T) is 
given by 


( 2 . 2 ) 


\/~D = [ao; 


a i, 


", 2ao] 


where do = 


y[D , the greatest integer < y[D\ the a,; 


(1 < i < r, where i and 


a, > 0, £ Z+) are the partial denominators of the continued fraction expansion; the 
finite sequence ai,...,a r is symmetric (so a,j = a r -j+i for 1 < j < |_r/2j , with 
j > 0, £ Z + ); and the bar indicates the period of length r + 1. We use then the 
algorithm given by Sierpinski [13] to calculate the continued fraction development 
of y/~D, that can be summarized as follows. Let bi and c, > 0, £ Z, Xi £ R. irrational 
numbers, dj = [%i\ for i > 1, the partial denominators ai of the continued fraction 
are then successively 


(2.3) d 0 


Vd = 


1 Ci 

do H-— do H—7=- 

xi VD + bi 


, bi = d 0 , Ci = D - b\ 


dj = L^iJ 


Xi = a i 


1 

-= di 

1 


VD + bi 



= d i 


Ci- 1-1 

VD + bi+i 


(2.4) 


bi -\-1 — did bi , Ci -\-1 — 


D ^ 6? +1 


(see [13], p.313, for more details). Therefore dj_i, bi and d can be calculated 
successively by 


(2.5) 


Oj-l = 


do + bi- 1 
d-l 


bi = 


ld — 1 bi — 


1 7 Ci - 


D-b'j 

d-l 


with bo = 0 and Co = 1, until dis found such as dj = d r +i = 2do- 
The fundamental solution (XpYi) of the simple Pell equation (1.2) is found for 
a particular value of D by computing the r th convergent ( p r /q r ) of the continued 
fraction [do; di,..., a r , 2do] of \[D. The convergents pt and q, can be found by the 
recurrence relations 


(2.6) pi = ciiPi-i + pi-2 ; ® = CLiqi-i + qi - 2 

with p -1 = 1, po = do, q ~i = 0, qo = 1. The fundamental solution (Xi, Yi) is then 
( p r ,q r ) if r = l(mo<i2), or (p 2 r+i, 92 r+i) if r = 0(mod2), with X\ and Yi forming 



ON SOME GENERAL SOLUTIONS OF THE SIMPLE PELL EQUATION 


3 


an irreducible fraction (X 1 /Y 1 ), where, to recall, r + 1 is the period length of the 
continued fraction expansion of yfl 5. 


3. Theorems on general solutions of the simple Pell equation 

A first theorem is demonstrated allowing to find the convergents ( p r /q r ) of \J~D for 
D = (/ (u)) ±2 Q n where n > 0 and a > 0, £ Z, and f (u) is any polynomial function 
of u £ Z, taking only positive integer values and such that / (u) = 0 ( mod (2“' ln )) 
for Vu £ Z. 


Theorem 1. For squarefree D, n and r > 0, £ Z and a > 0, £ Z, let f (u) be any 
polynomial function of u £ Z such that f (u) > 0, £ Z for \/u £ Z and f is written 
instead of f ( u ) /or convenience; if 

(3.1) D = f±2 a n 
and 

(3.2) / = 0 (mod (2“ -1 n)) 

1/ien the convergents ( p r /q r ) of \fT) are 

(3.3) ^ 

»• V 2=£; 

and for the case of the minus sign in front of 2 a n in (3.1), f 2 > 2“n and a and n 
not taking simultaneously the values a = 0 and n = 1. 

This theorem is demonstrated by calculating the continued fraction development of 
\fD following Sierpinski’s algorithm (2.3) to (2.5) and the convergents ( p r /q r ) of 
\fT) by (2.6). 


Proof. For squarefree D , n and r > 0, £ Z; a > 0, £ Z; / («) any polynomial of 
u £ Z taking values f = f (u) > 0, £ Z. 

(i) For the case of the plus sign in front of 2 a n in (3.1), i.e. 

(3.4) D = f 2 + 2 a n 
one has 

(3.5) f 2 < f 2 + 2“n < (/ + l) 2 
under the condition that 


(3.6) 

It yields successively a 0 
the condition, if a > 0, 

(3.7) 

or if a = 0, 


/ > 2“ _1 n 


\Jf 2 + 2 a n =/,&i = /, ci = 2“n; ar 


/ 


with 


/ = 0 (mod (2“ 1 n)) 


(3.8) 

It yields then & 2 

(3.9) 


2/ = 0 (mod n) 

/, c 2 = 1; a 2 = 2/ = 2ao, giving the continued fraction 

/ 


// 2 + 2“n = 


/; 


la-1, 


, 2 /, 
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with r = 1 , which, with ( 2 . 6 ), yields immediately 

P L= ( A + 1 " 

<27 \ 


(3.10) 




(ii) For the case of the minus sign in front of 2 a n in (3.1), i.e. 

(3.11) D = f 2 — 2“n 

and a and n not simultaneously a = 0 and n = 1 , one has 

(3.12) (/ - l ) 2 <f- 2“n < f 
under the condition that 

(3.13) / > 2“ -1 n 


v.f 2 - 2“nJ =(/-!), h = (/-!), 


It yields successively ao = 

ci = ( 2 / — 2 Q n — 1 ); ai = 1 with the condition 
(3.14) / > 2“ra 

One has then 62 = (/ — 2“n), C 2 = 2 Q n; 02 = ~ w ith the condition (3.7) 

if a > 0 or (3.8) if a = 0. It yields then 63 = (/ — 2 a n), C 3 = (2/ — 2 a n — 1); 
«3 = 1 , 64 = (/ — 1 ), C 4 = 1 ; <24 = 2 (/ — 1 ) = 2 ao, giving the continued fraction 

(/- 1 );!, (^-2).l,2(7-1), 


(3.15) yfp - 2“n = 
with r = 3, which, with (2.6), yields immediately 

(3.16) 


Pr 

<27 


) _ 1 
2° - 1 n ± 


/ 

2 “- 1 » 


□ 


Note that the conditions (3.2) and / > 0 means that 3 k > 0, G Z such that 
/ = 2 a ~ 1 kn, yielding simpler expressions of D and of the convergents ( p r /q r ) of 
\fD as 

(3.17) D = 2“n (2 a “ 2 fc 2 n ± l) 

(3.18) jV 7 2 °-‘Un±l \ 

q r \ k ) 

The following corollary give particular relations deduced from this Theorem, some 
of them being already given by other authors (see e.g. [13]). 


Corollary 2. For squarefree D, d, m, n and r > 0, £ Z; a > 0, £ Z; and for appro¬ 
priate conditions as in Theorem 1, the following relations hold for the convergents 
(p r /q r ) of \/D : 

(i) for D ={d 2 - l), 


(3.19) 


Pr_ 

Qr 


(ii) for D = {d 2 + l), 
(3.20) 


Pr _ { 2 d 2 + 1 \ 
<27 V 2d ) 
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Pr_ 

Qr 


Pr_ 

Qr 


Vr_ 

Qr 


d 2 ± 1 


t- ± 1 


2 a ~ 1 n 


(in) for D = (d 2 ± 2), 

(3.21) 

(iv) forD = m 2 ( d 2 ± 2 a n), 

(3.22) 

(v) for D = md /3 ) 2 ± 2 a n 

(3.23) 

(vi) forD = n ( nd 2 ± l), 

(3.24) 

(vii) forD = d ( m 2 d ± 2), 

(3.25) 

(viii) for D = d (d + 4), d = 0 ( mod 2) and d>2, 


d 

2 a-1 mn 


± 1 

-* 1 n 

mil- 1, 

2 a ~ 1 n 


Pr_ _ | 

( 2nd 2 ± 1 

q r 

V 2d 

Pr _ 

/ ?n 2 d ± 1 

Qr 

\ m 


(3.26) 


Pr 

Qr 


(d+ 2) 2 


- 1 


d+2 

2 


(ke) for D = ^(d ± 2) 2 ±4^, d = 0 (mod 2) and d> 2, 


(3.27) 


Pr 

Qr 


(d± 2) 2 


± 1 


d±2 

2 


Proof. For squarefree H, d, d!, m, n, n' and r > 0, G Z; t, a and /3 > 0, G Z; let 
(3.28) / (d) = (md^ ± t ) 

be a polynomial of d such that / (d) > 0, G Z. One has then the following from 
Theorem 1, yielding the respective convergents by (2.6): 

(i) let t = a = 0, and m = n = /? = 1 in (3.28) and (3.1), yielding /(d) = d 
and D = d 2 — 1. Then, in the proof (ii) of Theorem 1, <22 = 2 (d — 1)) = 2ao for 
d> 2 , and the continued fraction reduces to d 2 — 1 = [(d — 1 ); 1 , 2 (d — 1 ),...], 
with r = 1, yielding directly (3.19) by (2.6) 2 . Note that (3.3) yields a following 
convergent for r’ = 2r + 1 = 3, which provides a next solution to the Pell equation. 

(ii) let t = a = 0, and m = n = /3=1 in (3.28) and (3.1), yielding / (d) = d and 
D — d 2 T 1; then (3.3) yields immediately (3.20); 

(iii) let t = 0, and m = n = a = /3 = lin (3.28) and (3.1) yielding /(d) = d 
and D = d 2 ± 2, giving immediately (3.21) from (3.3), with, for the minus sign, 


2 1n [13] p.321, it is erroneously stated that “Va 2 + 1 = [a; a, 2a,.. .] (instead of 
[a; 2a,...]) for any natural number a” and “ \J (na 1 ) 2 -f I (instead of \J ( na ) 2 + 1 ) = 
[(na — 1); 1 , 2n — 2, 1 , 2 (na — 1),... ] for natural numbers a and n > 1” 
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the condition that d > 3; if d = 2, the continued fraction (3.15) reduces to \/2 = 
[1; 2,...] with r = 1, yielding (pi/qi) = (3/2) by (2.6). 

(iv) let t = 0, /3 = 1 and n = n'm 2 in (3.28) and (3.1), yielding / (d) = md and 
D = m 2 ( d 2 ± 2 a n') and dropping the then (3.22) is immediate from (3.3) with 
the condition that d = 0 (mod ( 2 “ _ 1 nm)); 

(v) let t = 0 in (3.28), yielding f (d) = md^ 3 ; then (3.23) is immediate from (3.3), 
with the condition that mdP = 0 ( mod ( 2 “' lf 0 ); 

(vi) setting t = a = 0, /3 = 1 and m = n in (3.23) yields immediately (3.24); 

(vii) let t = 0 and j3 = 1 and either a = 1 and n = d, or a = 0 and n = 2d, then 

(3.25) is immediate from (3.3); 

(viii) and (ix) let m, n and /3 = 1, t and a = 2 , and d = 0 (mod 2) with d > 2; then 

(3.26) and (3.27) are immediate from (3.3). □ 


Another theorem is demonstrated allowing to calculate the convergents ( p r /q r ) of 
\f~D for D = (/ (u)) 2 ±2 a n where / (u) is any polynomial function of u, taking only 
positive integer values and such that f (u) = (2 a ~ 2 n) (mod (2 a- 1 n)) for Vu G Z, 
with n > 0 and a > 0, G Z. 


Theorem 3. For squarefree D, n and r > 0, G Z and a > 0, G Z, /et / (u) be any 
polynomial function of u G Z smc/i that f (u) > 0, G Z for Vu G Z and f is written 
instead of f ( u ) /or convenience; if 

(3.29) D = f 2 ± 2“n 

and 


(3.30) / = (2° 2 n) (mod (2 a 1 n)) 

t/ien t/ie convergents ( p r /q r ) of ypD are, if f = 0 (mod 2 ), 


(3.31) 

and, if f = 1 (mod 2 ), 

(3.32) 


Pr 

dr 


/ 2 (/ 2 ±2 a n) 

—-r L 

22a—3 n 2 


Pr 

dr 


/ 2 (/ 2 ±3x2°~ 2 n) 


2 3 “- 5 n 3 


± 1 


/(/ 2 ±2°~ 2 ra)(/ 2 ±3x2°~ 2 n) 


23a —5 n 3 


under the conditions for D = (/ 2 + 2 a n) and f = 0 or 1 (mod 2), 

(3.33) if 1 < n < 2 2-a t/ien / > 2“ _1 n 

(3.34) if n > 2 2 ~ a then f>( 3 x 2 a " 2 n - l) 

for D = (f 2 — 2 a n) and f = 0 (mod 2), 

(3.35) if a > 2, or if a = 0or land 3n > 2 3_a t/ien / > 2 a n 

(3.36) */ a = 0 or land 3n < 2 3_a t/ien / > (2 a_2 n + 2) 

for D = (/ 2 — 2 a n) and f = 1 (mod2), 

(3.37) */ n > (2 2-a + 1) then f > 3 x 2““ 2 n 

(3.38) */ 1 < n < 2 2_a then f > 2 a ~ 3 n + 3 
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Proof. For squarefree D, n and r > 0, G Z; a > 0, G Z; / (u) any polynomial of 
u G Z taking values f = f (u) > 0, £ Z, let 

(3.39) / = (2 a ~ 2 n) ( mod (2“ _ 1 n)) 

Obviously for a > 3, f = 0(mod2) always, while for 0 < a < 2, / > 0, G Z 
if n = 0 (mod2 2 ~ a ) and / = 0 (mod 2) if n = 0 (mod 2 3_Q ), or / = 1 (mod 2) if 
n = 2 2 ~ a (mod 2 3- “). 

(i) For the case of the plus sign in front of 2“?i in (3.29), i.e. 

(3.40) D=(f 2 + 2“n) 
one has 

(3.41) / 2 </ 2 + 2“n<(/ + l ) 2 
under the condition that 


(3.42) 


/ > 2““ 2 n 


It yields successively do 
&2 = (/ - 2 “ _ 1 ?r), c 2 = 

(3.43) 


= |y/ 2 + 2“nJ 

(/- 2 “- 2 n + l) 


= fi h = /, ci = 2“n; m = 
a 2 = 1 with the condition that 


/ > 2 a ~ 2 n + 1 


Then & 3 = (2“ 2 n + 1), c 3 = (/ + 2° 2 n - 1); a 3 = 1, 64 = (/ - 2), c 4 = 4; 
/“I 

a 4 = 

(i.l) If / = 0 (mod 2 ), a 4 = (^), = (/ - 2), c 5 = (/ + 2 “~ 2 ?z - l); a 5 = 1 

with the condition (3.43). Then b 6 = (2 a ~ 2 n+ l), cq = (/ — 2“ _2 n + l); tig = 1 
with the condition that 


(3.44) 


/ > 3 x 2 a ~ 2 n — 1 


and 67 = (/ - 2 “ 1 n), c 7 = 2 “n; a 7 = = /, c 8 = 1 ; a 8 = 2 / = 2 a 0 , 

yielding the continued fraction 
(3.45) 

v7 2 + 2 “n = 


/; 


/ - 2 a ~ 2 n 


2 a ~ 1 n 

with r = 7, which by (2.6) yields 


, 1 , 1 , 


/ — 2 


, 1 , 1 , 


/ - 2 Q “ 2 n 
2“" 1 n 


, 2 /,... 


(3.46) 


P7 

97 


'/ 2 (/ 2 + 2 “») , 1 

—-•" 1 

/(/ 2 +2°~ 1 nT 

22a —3 n 2 


(i.2) If / = 1 (mod 2), a 4 = Then, b 5 = f, c 5 = 2 a 2 n ; a 5 = (5 

&6 = /, c 6 = 4; a 6 = (^); h = (f - 2), c 7 = (/ + 2 a ^ 2 n-l); a 7 = 1 with 
the condition (3.43), &g = (2 a-2 n + l), eg = (/ —2 a_2 n + l); as = 1 with the 
condition (3.44), 69 = (/ - 2“ _ 1 n), c 9 = 2“n; a 9 = 6 10 = /, c 40 = 1; 
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aio = 2/ = 2ao, yielding the continued fraction 

'f — 2 a ~ 2 n 


\// 2 + 2 a n = 


/; 


2 a ~ 1 n 


1 , 1 , 


/-I 


/ 


’ V 2 a ~ 3 n J ’ 


(3.47) 


/-I 


1 i | f ~ 2 °~ 2 ” 
’ ’ ’ 1 2“- 1 n 


, 2 /,.. 


with r = 9, which by (2.6) yields 


(3.48) 


P9 

99 


/ 2 (/ 2 +3> 


23a — 5 n 3 


+ 1 


/(/ 2 +2°~ 2 ri.)(/ 2 +3x2°- 2 n) 

23a — 5 n 3 


Summarizing the conditions (3.42) to (3.44) for / = 0(mod2 ) or / = l(mod2), 
one has that 

- if 1 < n < 2 2 ~ a , then 

(3.49) / > 2 Q_1 n 

- if n > 2 2_a , then 

(3.50) / > 3 x 2 a ~ 2 n — 1 

(ii) For the case of the minus sign in front of 2“n in (3.29), i.e. 

(3.51) D = (/ 2 - 2 Q n) 
one has 

(3.52) (/ - l ) 2 < / 2 - 2“n < / 2 

giving the same condition as (3.42). It yields successively ao = \Jf 1 — 2 Q nJ = 
(/ — 1 ), 6 i = (/ — 1 ), Ci = ( 2 / — 2 Q n — 1 ); ai = 1 with the condition that 

(3.53) / > 2 a n 

Then, b 2 = (/ - 2“n), c 2 = 2 Q n; a 2 = w ), 63 = (/-2 a_ 1 n), c 3 = 

(j _ 2 a_2 n — l); a 3 = 2 with the condition that 

(3.54) / > 2 a ~ 2 n + 2 

and 64 = (/ - 2 ), c 4 = 4; a 4 = • 

(ii.l) If / = 0 (mod 2), a 4 = Then, b 5 = (/ - 2), c 5 = (/-2 a_ 2 n-l); 

as = 2 with the condition that 

(3.55) / > 3 x 2 a ~ 2 n 

Then, & 6 = (/ — 2“ -1 n), c 6 = 2“n; a 6 = (^1^), 6 7 = (/- 2“n), c 7 = 
(2/ - 2“n - 1); a 7 = 1, 6 8 = (/ - 1), c 8 = 1; a 8 = 2 (/ - 1) = 2a 0 , yielding the 
continued fraction 


v7 2 -2 


— 


(/- i);i> 


/ - 3 x 2 a ~ 2 n 


2 a ~ 1 n 


/ — 2 


(3.56) 


2, 


/ - 3 x 2“" 2 n 
2“- 1 ?r 


, 1,2 (/ - 1 ), 
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with r = 7, which by (2.6) yields 


(3.57) 


P7_ 

17 


22a —3 n 2 I -L 

/(/2-2°-i^T 

2 2a —3 n 2 


Summarizing the conditions (3.42) and (3.53) to (3.55) for / = 0 (mod 2), one has 
that 

- if a > 2, or if a = 0 or 1 and 3n > 2 3_a , then 

(3.58) / > 2 a n 

- if a = 0 or 1 and 3n < 2 3_a , then 

(3.59) / > 2 a ~ 2 n + 2 

(ii.2) If / = 1 (mod2), a 4 = h = (f ~ 4), c 5 = (2/ - 2 a ~ 2 n - 4); a 5 = 1 

with the condition that 

(3.60) / > 2 a_1 n + 1 

and b 6 = (f-2 a ~ 2 n), c 6 = 2 a ~ 2 n; a 6 = h - 

(2/ + 2 a ~ 2 n — 4); a 7 = 1 with the condition that 

(3.61) / > 2 a ~ 3 n + 3 

Then & 8 = (/ - 4), c 8 = 4; a 8 = b 9 = (f - 2), c 9 = (/-2“ _2 n-l); 

ag = 2 with the condition (3.55), &io = (/ — 2 Q_1 n), cio = 2“n; 
aio = = (/ - 2“^)- cn = (2/ - 2 Q n - 1); o n = 1, b 12 = (f - 1), 

C12 = 1; ai2 = 2 (/ — 1) = 2ao, yielding the continued fraction 


( V-Sn" ); = (/ — 2“ 2 n), c 7 = 


\// 2 - 2“n = 


(3.62) 


(/-i);i» 

/ - 2 a ~ 2 n 
2 a ~ 3 n 
/ - 3 x 2 a ~ 2 n 


f-3x2«~ 2 n\ (f- 3 

> 


1, 


2“ _1 n 
/ — 3 


, 1 , 


2“ _1 n 

with r = 11, which by (2.6) yields 

f 2 (f 2 -3x2°‘~ 2 n) 


, 2 , 


) 1 ) 2 (/ — !),••• 


(3.63) 


Pn 

9n 


23q-5„3 


^ 1 


/(f —2° — 2 n)(/ 2 —3x2“ — 2 f 
— 5 n 3 


Summarizing the conditions (3.42), (3.60), (3.61) and (3.55) for / = 1 (mod 2) and 
as 0 < a < 2, one has that 

- if n > 2 2_a + 1, then 

(3.64) / > 3 x 2 a ~ 2 n 

- if 1 < n < 2 2 ~ a , then 

(3.65) / > 2 a ~ 3 n + 3 


□ 
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Like above, note that the conditions (3.30) and / > 0 means that 3 k > 0, € Z such 
that / = 2“ -2 (2k + 1) n, yielding from (3.29) 

(3.66) D = 2“n (V" 4 (2fc + l) 2 n ± l) 


and for the plus sign in front of 1 in (3.66), conditions (3.33) and (3.34) reduce to 
k> 1, and for the minus sign in front of 1 in (3.66), for / = 0 (mod2), (3.35) and 
(3.36) reduce respectively to k > 0 and k > 2, and for / = 1 (mod2), (3.37) and 
(3.38) reduce respectively to fc > 1 and k > 2. Expressions simpler than (3.31) and 
(3.32) for the convergents (p r /q r ) of y/D are then 


(3.67) 


Pr_ 

q r 


V" 1 (2k + l) 2 n (V “ 4 (2k + l) 2 n ± l) + l' 
(2k + 1) ( 2 a ~ 3 (2k + l) 2 n ± l) 


if f = 2“ 2 (2k + l)ns0 (mod 2) (i.e. for either a > 3, or if a = 0, n = 0 (mod 8), 
or if a = 1, n = 0 (mod A), or if a = 2, n = 0 (mod 2)), and 


(3.68) — = 

Qr 


2 a ~ 3 (2k + l) 2 n (2 a ~ 2 (2k + l ) 2 n ± 3) ± 1 
i (2k + 1) (2“~ 2 (2k + l) 2 n ± l) (2“" 2 (2k + if n ± 3 ) 


if / = 2“ -2 (2k + 1) n = 1 (mod 2) (i.e. for a < 2, if a = 0, n = 4 (mod 8), or if 
a = 1, n = 2 (modi), or if a = 2, n = 1 (mod 2)). 

Several particular relations can also be deduced from the general cases given in this 
Theorem. 


Corollary 4. For squarefree D, d, m, n and r > 0, £ Z, d = 1 (mod2); a > 0, G Z; 
and for appropriate conditions as in Theorem 3, the following relations hold for the 
convergents (p r /q r ) of y/T) : 

(i) for D = 4d (d ± 2) and, in the case of the minus sign, d > 3 , 

( 3 . 69 , *-< 2 ^* 2 ) + !' 


<27 


d ± 1 


(ii) for D = 4 (d 2 ± 2) and, in the case of the minus sign, d > 3, 

p r __ / 2d 2 (d 2 ± 2) + 1 \ 


(3.70) 


Qr 


d(d 2 ± 1) 


(iiij for D = 16n (nd 2 ± l) and, in the case of the minus sign, nd > 1, 

p r f 8nd 2 (nd 2 ± l) + l\ 


(3.71) 




d (2nd 2 ± 1) 


(iv) for D = An (nd 2 ± 2) and, in tde case of the minus sign, nd > 2, 

p r f 2nd 2 (nd 2 ±2) + 1 \ 


(3.72) 


Qr 


d (nd 2 ± 1) 


(v) for D = (d 2 ± 4) and, in the case of the minus sign, d > 3 

' d 2 (d 2 ±3) 2 

(3.73) 


Pr_ 

Qr 


±1 

d(d 2 ±l)(d 2 ±3) 

2 
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(vi) for D = n (nd 2 ± 4) , n = 1 (mod 2) and, in the case of the minus sign, nd 2 > 2, 


(3.74) 


Qr 


/ nd 2 (nd 2 ±3) 2 \ 

I <1 ( n I IS ( n d 2 :1 3) \ 


Proof. For squarefree D, d, m, n, n' and r > 0, £ Z, d = 1 (mod 2), a and /3 > 0, £ 
Z; let 

(3.75) / (d) = md® 

such that / (d) > 0, £ Z. One has then the following from Theorem 3, yielding the 
respective convergents by (2.6): 

(i) immediate for m = 2, (3 = 1, a = 3 and n = d in (3.75) and (3.31). 

(ii) immediate for m = 2, (3 = 1 ,a = 3 and n = 1 in (3.75) and (3.31). 

(iii) immediate for m = 2 n, (3 = l,a = 3 and n = ‘In' and dropping the prime sign 
in (3.75) and (3.31). 

(iv) immediate for m = 2 n, (3 = 1 and a = 3 in (3.75) and (3.31). 

(v) immediate for m = 1, /3 = 1, a = 2 and n = 1 in (3.75) and (3.32). 

(vi) immediate for m = n, (3 = 1 and a = 2 in (3.75) and (3.32). □ 


For the cases (i) to (iv), it is easy to see that D is always 4 (mod8). From case 
(iv), generalizing cases (i) to (iii), D = 4n (nd 2 ± 2) = ^(2?rd) 2 ± 8 iij and it is easy 
to see that, as d = 1 (mod 2), if n = 0 or 1 (mod 2) then (2nd) = 0 or 2 (mod 4), 
i.e. solutions will be found for even integer squares (2nd) different from powers of 
2 ((2nd) ^ 2“), plus or minus even or odd multiples of 8 depending on the 0 or 
2 (mod4) value of 2nd. For example, for n = 1 (mod2), D = 6 2 ±8, 10 2 ±8, 14 2 ±8, 
... and also D = 30 2 ± 8, 30 2 ± 24, 30 2 ± 40, ...; for n = 0 (mod2), D = 12 2 ± 16, 
20 2 ± 16, 24 2 ± 32, 28 2 ± 16, ... 

The case (v) of D = (d 2 + 4) with d = 1 (mod 2) is interesting as it gives the solu¬ 
tions for all the values of D of the form D = (8A (k) + 5) = ((2k + l) 2 + 4^j with 

k > 0, £ Z, i.e. D = 5 (mod8) and = A (A:), where A (k) is the triangular 

number of k. Similarly, the case (v) of D = (d 2 — 4) with d = 1 (mod 2) gives the 
solutions for all the values of D of the form D = (8A (k) — 3) = ((2k + l) 2 — 4^j 

with k > 1, £ Z, i.e. D = —3 (mod 8) and = A ( k ). For example, D = 3 2 +4, 

5 2 ± 4, 7 2 ± 4, 9 2 ± 4, ... 

The case (vi) is also interesting as it generalizes the previous case (v). D can in¬ 
distinctly be written as D = n (nd 2 ± 4) = (8A ( wd 2 ~ 1 ) + 1 ± 4rc), which is always 
D = 5 (mod 8) as d and n = 1 (mod 2). For example, D = 9 2 ±12, 15 2 ±12, 15 2 ±20, 
21 2 ± 12, 21 2 ± 28, ... 


4. Conclusions 

The case considered in these two theorems of any polynomial function / (u) of u £ Z 
can easily be extended to any function g (x) of x £ R or h (z) of z £ C as long as 
g (x) or h (z) > 0, £ Z for Vx £ R or Mz £ C. 

These two general theorems and the various relations in the two corollaries allow 
the immediate calculation of the fundamental solutions of the simple Pell equation 
for most of the values of D. A quick survey shows that these two theorems provide 
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solutions for all squarefree integer values of D for D < 10, for 72.2% of the values 
of D for D < 10 2 , for 35.4% for D < 10 3 , and for 15.7% for D < 10 4 . 
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